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One dimensional nonlinear plasma wake-fields excited by a single bunch and by series 
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density are investigated. The obtained results shows dependence of excited potential 
electric fields on bunches parameters and allows to choose these parameters optimal. 
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1 Introduction 



The electromagnetic waves, excited in plasma by charged bunches, can be used both for 
focusing of bunches and for charge acceleration [I]. In the relativistic bunch case an ampli- 
tude of excited one dimensional nonlinear plasma waves can essentially exceed the conven- 
tional wavebreaking field E* = m e VQUJ p /e (t>o is the phase velocity, uj p = (ATcn e 2 /m e ) 1 ^ 2 
is the electron plasma frequency, Uq- plasma electron density in equilibrium) and reaches 
of the value E = [2(7 — l)] 1 ^ 2 , where 7 = (1 — /3 2 ) _1,/2 , (3 = Vq/c, and E is normalized on 
E*. Acceleration rate in the wake- fields can reaches of a few GeV/m, that is much more, 
than was reached on conventional accelerators. 

The one dimensional nonlinear wake-fields theory was developed in series of works 
(see e.g. ||-[f§ and references therein). Maximum of the accelerated electrons energy, as 
it was shown in ||, can reach of value 4m e c 2 7 3 . 

In this paper bunches with the different parameters investigated analytically and by 
numerical simulation (including nonuniform and positively charged bunches). 



2 One Dimensional Nonlinear Wake-fields 

Consider a cold electron plasma with moving through the plasma, in Z direction charged 
bunch (or train of bunches), which transverse sizes one can consider infinite. The sub- 
ject of our investigation is the steady wakefields excited by such bunches. Plasma ions 
assumed immobile because of their large mass (m e /mj < 1). As usually we also as- 
sume, that field in front of the bunch (or bunches) is absent, potential and strength of 
the electric field are continuous in all space. Then, from the continuety equation, rela- 
tivistic equation of motion for plasma electrons and the Poisson equation one can obtain 
following equation for steady fields ||,|§ 

g + ^V(l- (f2 _^_ 2)1/ j +/? 2 a(,)=0, (1) 

where F = 1 + \e\ip/m e c 2 > 7" 1 , if is the electric potential, z = k p (Z — v t), k p = uj p /v 
(here the phase velocity vo is equal to bunch velocity), a(z) = (q/\e\)nb(z)/n , e is the 
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electron charge, q is the charge of the bunch particles, n\,{z) is their density. Normalized 
on E* the strength of electric field obeys the formula E z (z) = E — —{l/f3 2 )dF/dz. The 
plasma electrons dimensionless velocity as function of F is 

fie = \P — (F 2 — 7~ 2 )^ 2 ]/ (fl 2 + F 2 ). (2) 
Inside uniform bunches and outside of bunches Eq. ([!]) can be rewritten in the form 

d 2 F dU_ Q 

U = /? 2 [( 7 2 + a)F - /? 7 2 (F 2 - 7 ~ 2 ) 1/2 ]. (3) 

Formally Eq. (|3|) describes one-dimensional motion of a particle in a field with potential 
U(F). Analysis of the function U(F) allows to determine both qualitative behaviour of 
the field and some quantitative values (such as electric field amplitude in the bunch, the 
wave-breaking field, maximum value of electric potential and other). Moreover, it can 
be obtained analytical solution of Eq. (|]), which includes the elliptic functions [|],|5|-[0]. 
However, even in the case of uniform bunch, the analytical description of the field behind 
bunch in general case is practically impossible. For nonuniform bunches the problem 
analytically is not solved yet. 

Outside the bunch (where at(z) = 0) "potential" U(F) has one minimum, therefore 
the field is periodical and for fixed 7 fully determines by electric field amplitude E mp < E. 
It is necessary to note that Eq. (|]) with a = describes also the wake wave, excited by 
other way (for example, by laser pulse). 

Nonlinear wake wavelength A p increases as the amplitude E mp increases. When 7 ^> 1 
the dependence A p (E mp ) practically not depend on 7 and is presented on Figure 1 (notice, 
that according to accepted in ([!]) variables, value A p = 2tt corresponds to the linear 
approach). Near by the "breaking" (E mp w E) the wake wavelength in the case of 7 ^> 1 
is approximately equal to 4(27) x / 2 |J, || and maximum value of the dimensionless electric 
potential is F mp m (l+/5 2 )7 ||. In general case F mp increases as E mp increases (see Fig.l). 

The energy gain of electrons (or positrons) accelerating in the wake wave can reach 
of value 4m e c 2 7 3 ||. Really, the relativistic equation of motion of accelerating electron 
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in the frame of reference moving with the wave velocity is 

d{P' a i a )/dt' = cdF'/dZ'. (4) 

Taking into account, that dZ' = cf3' a dt' and f3' a d{f3' a / y' a ) = dj' a from @) follows 7^ = 
7a(0) + F' — -^'(0), where 7^(0) and -^'(0) are initial values. In the laboratory frame 
of reference 7 a one can obtain using known relativistic transformations tp = jip' and 
7a = 77 ^(l + /3/3'J. In the case 7 > 1 and (3 a (0) > (3, choosing F(0) « I/7 and 
F ~ F mp fa 2 7 (that corresponds to E mp « E) obtain 7a « 4 7 3 . In this case acceleration 
occur on the half of the wavelength;the acceleration length is proportional to 7 5 / 2 ||. For 
example, when n = 3 ■ 10 12 cm~ 3 and 7 = 10, the acceleration length is Z a ~ lm, A p « 
5cm and (•y a )max ~ 4000. The result (7 a )max ~ 7 3 is consequence both of correlations 
E m p ~ 7 1//2 ,A p ~ 7 1 / 2 and relativistically large mass of the accelerating particle (i.e. 
/3 a « /3 « 1), that causes long keeping back of the particle in accelerating phase. For 
arbitrary wake wave amplitude the maximum value F mp can be obtained from Fig.l. 

For the large 7 the acceleration length can surpasses the laboratory plasma sizes. 
In the cosmos (for example, in the blanket of supernova), where the plasma can be 
considered unlimited, acceleration in the strong plasma fields described here can originate 
high energy cosmic rays. 

3 Single Bunch 

Consider a single uniform bunch :nb = const, when — d < z < and = elsewhere. 
The field inside the uniform bunch describes by Eq. (^). "Potential" U(F) is maximum 
when F = I/7 and is minimum when 

F = F c = | 7 2 + a|/( 7 r 1/2 ), T = 7 2 (1 + 2a) + a 2 . (5) 

From ([5]) follows, that in the case —1/(1 + /?) = a.\ < a < 00 function U(F) has minimum 
and the field inside bunch is periodical;/ 71 changes in the limits 

l/7<^< [(7 2 + «) 2 + /?V]/7r. (6) 
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When a = ati,U(F) monotonously decreases as F increases and tends to zero when 
F — > oo. In this case the electric field strength inside the bunch grows from zero at the 
bunch head (z = 0) and tends to some constant value (see ([/])). When a < ai,U(F) 
decreases unlimitedly, E grows from the bunch head to its tail (z = —d), and for bunch 
length d more than some value d max the wake behind the bunch "breaks". 

When U(F(z = 0) = 1) > U rnax = U(F = I/7) (that takes place, when a > 7), 
in enough long bunch F becomes less than I/7, that is using in this paper cold fluid 
approach becomes inapplicable. Thus, in the cases —00 < a < a\ and a > 7 the uniform 
bunch length can not exceed some value d max ;the wake amplitude in the case of d = d max 
is equal to E. Dependence of the d max on a for 7 = 10 presented on Figure 2. 

In the case of ol\ < a < 7 the uniform bunch length can be arbitrary. The field inside 
the bunch is periodical and has amplitude ||, || 

E mb = [2{l+ a -T l ' 2 h)] 1 ' 2 /(3. (7) 

Dependence of the wavelength inside the bunch A& on a for 7 ^> 1 shown on Figure 3. 

In the linear approach E mp ~ a (\a\ <C 1) and differences between the electron 
bunch case and the positron one reduce to changing of electric field sign (saying electron 
or positron bunch we mean negatively or positively charged bunch, respectively). In 
nonlinear theory these differences are essential. The wavelength A& in the electron bunch 
increases as \a\ increases when a.\ < a < (see Fig. 3), and in the case 7 ^> 1 and a ~ 
—0.5 is nearly equal to 87 @, 0, |J, that is much more than linear plasma wavelength. 
In the positron bunch A^ decreases as a increases. Note, that difference of the nonlinear 
wavelength from linear one, in the case 7 ^> 1 is palpable even for |a| <^ 1. From (Q), 
(ID and the boundary conditions one can see, that F(z) > 1 in the electron bunch, and 
F(z) < 1 in the positron bunch. That is, according to (0), f3 e (z) < in the electron 
bunch and /3 e (z) > in the positron one. 

Equation(|l]) was investigated numerically for uniform, parabolic and linear charge 
density distributions in the bunch 



5 



a p = a p0 [l-(z-d/2) 2 ]/(d/2) 2 , 
ai = aio\z\/d, 
-d<z<0. (8) 

Bunches with profiles (||), but with the same length and same total charge are compared. 
That means |a p0 | — (3/2)|a M0 |> \ a io\ — 2|a n0 |- 

Typical dependence of the wake wave amplitude E mp on bunch length for distributions 
(H), in the case of periodical field inside the uniform bunch (i.e. when a± < a u o < 7) 
shown on Figure 4. For the uniform bunch this dependence is periodical (on Fig. 4 shown 
one period only). The wake wave amplitude behind the uniform bunch is maximum and 
equal to E mp = 2|a. lt o|7/r 1 / 2 || when d = (n + 1/2)A 6 (n is integer) and equal to 
zero when d = nh^. For the parabolic and linear profiles choice d « (3/4) Aj is optimal. 
One can see also, that in the case d > A& the "parabolic" bunches becomes less effective 
and the wake wave amplitude behind the bunch with linear profile weakly depends on d. 
For a short bunch (d < A&/2) E mp (d) almost not depends on the bunch profile. 

4 Trains of Bunches 

The results presented above demonstrates, that for the strong wake wave (with E mp ~ E) 
excitation it is necessary relativistic either negatively charged bunch with a < —0.5 or 
positively charged bunch with a > 7. However, on experiments density of bunches is 
often much less of the plasma density, i.e. \a\ 1. Therefore, it is naturally to try 
excite the strong waves by series of bunches. For the case of uniform bunches with 
the same densities this problem was investigated analytically in 0. The wake wave 
amplitude behind the train of uniform bunches is maximum when di = (rii + 1/2) A^ and 
lj = (rij + 1/2) A P j, where di is the length of i-th bunch, lj is the spacing between j-th 
and (j + l)-th bunches, A&j is the wavelength in i-th bunch, A p j is the wavelength behind 
j-th bunch, and rij are integer. The wake wave amplitude behind such optimized train 
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depending on number of bunch N is : 

E mp (N) « (2/p)8h(\a\/3N), \a\ < 1. (9) 

From (|^) follows, that the wave amplitude quickly grows as N increases and reaches the 
value E when 

N « (M/?)" 1 ln[( 7 + l) 1 / 2 ^ + ( 7 - l) 1 / 2 ^]. (10) 

The electric field in the optimised train, obtained by numerical simulation of equation (fj) 
presented on Figure 5. The numerical results showed high exactness of the expressions 
© and (0). 

Often in experiments using trains of bunches with d$ = d = const, lj — I — const and 
<C 1. Analytical discription of the fields excited in such trains is very difficult. 
The numerical simulation showed, that for profiles (HD in the case of d and I optimal for 
first bunch (see Sec. 3) the wake wave amplitude behind the train grows almost linearly 
up to N = iV* ~ 0.b/\a u o\.E mp (N*) increases as |a u0 | and (or) 7 increases. When the 
wave becomes essentially nonlinear (when N > N*) coherence of the waves exciting by the 
bunches breaking because of increasing of the wave length and the wave amplitude behind 
train decreases almost up to zero as N increases. Dependence E mp (N) for N 3> iV* is 
nearly periodical. 

It was ascertained, that trains of short bunches (d A&/2 and / ~ A p /2) are not 
effective as in this case E mp (N) < E mp {l). When d A b and I A p the train behaves 
as one bunch with averaged along the train length charge density. 

5 Conclusions 

The analytical and numerical results presented above throws light upon a number of 
questions of the theory of strong plasma waves excited by charged bunches. This results 
allow to choose parameters of the bunches optimal on future experiments. 

The peculiarity of wake waves excitation by positively charged bunches are revealed. 
In particular, it was discovered, that the nonlinear wavelength in uniform positively 
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charged bunch decreases as charge density increases. 

For the excitation of wake wave with E mp ~ E more suitable is a bunch, consisting 
of heavy particles (for example, protons), as light particles (electrons, positrons) lose 
considerable part of their energy on length comparable with the bunch length. This is 
easy to show by comparing of the bunch energy with the energy of wake wave. 

Note, that the results presented in this paper are also suitable for the case of two- 
dimensional cylindrical bunches with radius a when (■yc/upO) 2 <C 1 and r 2 Ca 2 ||. 
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Figure 1. Dependence of the wake wavelength and maximum value of the electric 
potential on wake wave amplitude E mp (in the dimensionless units); 7> 1 
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Figure 2. The maximum permissible bunch length d max depending on the bunch 
charge density (7 = 10). 
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Figure 3. Dependence of the wavelength in uniform bunch on the charge density. 
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Figure 4. The wake wave amplitude behind bunch depending on the bunch length 
(a u o = 0.1,7 = 10). 1-uniform bunch, 2-parabolic profile of the bunch charge density, 
3-linear profile case. 
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Figure 5. The electric field in optimized train of uniform bunches. 

a) The train of negatively charged bunches (a = —0.2,7 = 10). Inside bunches E > 0. 
Negative E corresponds to space between bunches. 

b) The train of positively charged bunches (a = 0.2,7 = 10). Inside bunches E < 
0. Positive E corresponds to space between bunches. 
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